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Abstract 

Analyzing T(nS) decays acquired with the CLEO detector operating at the CESR e + e~ collider, 
we measure for the first time the product branching fractions B[T(nS) — > jXb((n—l)Pj)] xB[xb(n— 
l)Pj) — > Xi] for n = 2 and 3, where X{ denotes, for each i, one of the fourteen exclusive light- 
hadron final states for which we observe significant signals in both Xb(^P.j) and Xb(2Pj) decays. 
We also determine upper limits for the electric dipole (El) transitions T(3S) — > 7Xb(l-Pj). 
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In the 31 years since the first observation of bottomonium we have learned a great deal 
about decays of the T(1S,2S,3S) resonances and transitions among them. Less is known 
about P-wave states because they are not produced directly in e + e~ collisions. The spin- 
triplet Xb mesons are produced copiously in electric dipole (El) transitions [lj, permitting 
the recent first observations of inclusive decays of Xb{ n Pj) to p(p) [1] and to open charm 0]. 
Nothing else is known about XbinPj) decays to non-66 states. Such processes are of interest 
both intrinsically and as clues in searching for states of mass ~ 10 GeV/c 2 via their exclusive 
decays. 

In this article we report the first observations of decays of XbO-Pj) and Xfe(2Pj) into 
specific final states of light hadrons, where the Xb(nPj) states are produced via T(2S) — ► 
7Xb(lPj) and T(3S) — > 7Xb(2Pj). We also determine upper limits on rates for the suppressed 
El transitions T(35) — > 7x b (lPj). 

We use the same T(nS) on-resonance data as in the analysis of Ref. Q, corresponding 
to N T(n5) = (20.82 ± 0.37, 9.32 ± 0.14, 5.88 ± 0.10) x 10 6 resonance decays for n = 1, 2, and 
3, respectively, collected by the CLEO III detector [Ej at the the Cornell Electron Storage 
Ring. Hadronic events are selected based on the criteria used in the analysis of Ref. 

Our signal events have the form T(nS) — *> 7X4, where Xi denotes a specific fully recon- 
structed final state. We allow a large variety of possibilities for Xi, but to keep the list finite 
and realistic we impose the following requirements. Each Xi consists of a combination of 
twelve or fewer "particles," where a "particle" is defined here to be a photon or a charged 
pion (tt ± ), kaon (K^), or proton (p/p). Each state Xi must have at least two charged 
"particles" and conserve overall charge, strangeness, and baryon number. We only consider 
modes in which photons other than that from the transition are paired into either 7r° or rj 
candidates, of which we only permit four or fewer. Neutral kaon decays into 7r + 7r~ are also 
considered. With these criteria, there are 659 separate final states, which act as the basis 
for our search. 

Photon candidates are taken from calorimeter showers that do not match the projected 
trajectory of any charged particle and which have a lateral shower profile consistent with 
that of an isolated electromagnetic shower. Each candidate for a 7r ± , K , and p/p must be 
positively identified as such by a combination of its specific ionization dE/dx, within 3a, 
where a refers to uncertainty due to measurement errors, and, when available, the response 
of a Ring Imaging Cherenkov system as in the analysis of Ref. (3)]. Candidates for 7r° and rj 
decays to two photons are allowed only if the photon pair mass is within 3a of the nominal 
7t° or rj mass. Kg — > 7r + 7r~ candidates, consisting of a pair of vertex-constrained oppositely 
charged tracks, are required to have effective mass within 3a of the nominal mass [I] an d to 
have a flight path before decay exceeding twice the longitudinal vertex resolution. 

We improve sample purity by constraining the transition photon plus the decay products 
of Xi to the initial T(nS) four-momentum with a 4C kinematic fit and requiring the fitted 
X 2 (4C)/dof < 5 as in Ref. [cjj]. The kinematic fit also allows us to improve the resolution on 
the invariant mass of the Xi by using fitted, instead of measured, four-momenta: we denote 
this mass by M inv . 

Figure [JJ shows fits to the M inv distribution of the sum of all 659 modes in (a) T(2S) 
and (b) T(35) data. The natural Xb{nPj) widths Q are expected to be much smaller than 
the resolution (~ 5 MeV) of the transition photon, which has a mostly Gaussian line shape 
and a low-energy tail induced by energy leakage out of the crystals used in the algorithm. 
This Crystal Ball line (CBL) shape is discussed in more detail in Ref. To fit the M- mv 
distribution in Fig. [TJ we use a "reversed" CBL shape with the asymmetric tail on the 
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FIG. 1: Mi nv spectra based on (a,c) T(2S) and (b,d) T(35) data for the sum of all 659 modes (a,b) 
and the 14 selected modes (c,d). The observed peaks are consistent with the transitions T(25) — > 
7X;,(lPj) and T(35) — > 7Xb(2Pj). Fitted backgrounds are represented by dashed histograms, 
fitted Xb{ n P.j) peaks are represented by dotted lines, and sums of fitted signals and background 
are denoted by solid curves. 



high side instead of the low side of thepeak. The fitted masses of XbO-Pj) an d Xfe(2-Pj) 
are consistent with the known masses For background shapes, we use Mi nv spectra 
obtained with the same analysis procedure but based on Y(IS') data, shifted by differences 
in center-of-mass energy while floating normalizations. This procedure appears to represent 
backgrounds reasonably. Using low-order polynomials instead of T(IS') data to represent 
backgrounds, we obtain consistent results. 

With signal shapes, including central values, fixed by a fit to the sum of the 659 modes, 
we fit the unconstrained photon energy spectra for each mode with CBL shapes. We use 
unconstrained photon energy spectra because calorimeter resolutions are independent of 
the final states in ^(lP, 2P) decays. We then determine significances from the fit to each 

mode as \j —2 \n(C wo / C w ) where C wo and C w are likelihoods from fits without and with 
an allowance for signal. We determine the significance from simultaneous fits to the three 
peaks instead of determining the significance of individual x&(1Pj, 2Pj) peaks. We identify 
14 modes giving at least 5<r significance from both Xb(l-Pj) and Xb(2Pj) decays. 

On the basis of GEANT-based [9J signal and various background Monte Carlo (MC) 
samples for the 14 identified modes, the required limit on x 2 (4C)/dof is varied from its 
initial value of 5 in order to optimize signal sensitivity while reducing backgrounds. The 
optimum value for the 14 modes is found to be x 2 (4C)/dof < 3, and is adopted as our 
nominal value. As some modes show further improvement in sensitivity for x 2 (4C)/dof < 2, 
we also explore the choices x 2 (4C)/dof < 2 and < 4 in our study of systematic uncertainties. 

Fig. [T] shows M inv distributions of (c) T(2S) and (d) T(3S f ) data based on the sum of 
the 14 modes with our nominal selection criteria. The fitted backgrounds are T(IS') data, 
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FIG. 2: Fits to invariant masses of individual decay modes. For each mode X, (a) refers to 
T(2S) -> JXb(lP), XbO-P) -> Xi, while (b) refers to T(35) -> 7X&(2P), Xb(2P) -» X<. In labels 
for X, iir^ = i^g. Dotted lines represent fitted constant backgrounds while dashed lines are fitted 
Xb(nPj) signals of each J state. 



shifted as in Figs. Q](a,b). The fitted Xb{nPj) masses again are consistent with the known 
values [lj]. With the restriction of x 2 (4C)/dof < 3, T(2S,3S) decays in the 14 modes lead 
to roughly 40% of the total observed events in the 659 modes. 

To measure B[T(nS) — > ■JXb(( n — 1)Pj)] x &[Xb(( n — 1)Pj) Xi], where Xi is each of 
the 14 modes, fits to signal Monte Carlo samples for signals produced through transitions 
of T(nS) — > 7Xfe(( n — 1)Pj) are performed to M inv spectra. Once signal shapes are fixed for 
each mode, we perform fits to data. We fix the central values of ^(lP, 2P) masses according 
to world averages [lj]. Fitted x 2 (4C)/dof distributions for each mode of X&(1P> 2P) and each 
J are found to behave as expected from signal MC samples. Resultant fits to M inv spectra 
are shown in Fig. [2j For all cases, we use a constant (flat) background shape and fitting 
ranges of 9800-9950 MeV/c 2 in T(2S) data and 10180-10300 MeV/c 2 in T(3S) data. 

The major source of systematic uncertainty is found to be the effects on signal efficiency 
of possible intermediate states. We study this in T(25) — > 7X&(1P) and apply the result to 
T(3S) — > 7x 6 (2P) as well. In all our signal Monte Carlo samples, Xb decays were generated 
according to phase space. To estimate the systematic uncertainty due to the presence 
of intermediate states, we consider — > 7T ± 7r°, p° — > 7r + 7r~, <p — > K + K~, uj — > 7r + 7r _ 7r , 
r) -> tt+tt-tt , K*(892)° -> K^/K^tt , and fsf*(892) ± -> ^vr /^^. We find deviations 
in the efficiencies based on these modes to be as large as 18%; hence to allow for possible 
neglected intermediate states, we assign ±20% systematic uncertainty due to MC modeling 
of Xb(nP) decays. 

Other sources of systematic errors were found to be small in comparison with the pos- 
sible presence of intermediate states. In roughly descending order of importance, they are: 
kinematic fitting (7-14%); photon, 7r°, and charged track reconstruction (4-10%); parti- 
cle identification and Kg efficiencies (4-10%); statistical uncertainty on signal MC samples 
(2-8%); numbers of T(nS) (2%); cross feeds among our 14 signal modes (1%); fit ranges; 
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background shapes; bin width; signal widths; peak positions of Xb', trigger simulation; and 
multiple candidates. Systematic errors are added in quadrature mode by mode. They fall 
within a range of 23-30% for all modes. 

Table [I] shows efficiencies, yields, and signal significances for each of the 14 modes of 
T(2, 35") — + 7Xfe(l-Pj, 2Pj). Table HT1 shows the measured product branching fractions 
B[T{nS) 7Xb((n-l)Pj)]xB[xb{(n-l)Pj) -> X t ], forn = 2 and 3, in units of KT 5 . For all 
transitions with significance less than 3a, we set upper limits at 90% confidence level (C.L.), 
also shown in this table. Table [TTT1 shows the measured values of B[xb((n — l)Pj) — > XJ ob- 
tained using the values of B[Y{nS) -> VCb((n-l)Pj)] = (3.8±0.4, 6.9±0.4, 7.15±0.35)% for 
n = 2 and (5.9±0.6, 12.6±1.2, 13.1±1.6)% for n = 3, and for J = 0, 1, and 2 respectively Q 
whose uncertainties are also included in the systematic errors. 

As expected for particles of mass ~ 10 GeV/c 2 , exclusive decays are distributed over 
many final states. The values of B[xb((n—l)Pj) — > Xj\ listed in Table [TTTI are typically a few 
parts in 10 4 , suggesting that the decay modes of these 10 GeV particles are distributed over 
more than a thousand different modes, of which we have investigated 659. Several points 
are worth noting. 

(1) The mode with the largest branching ratio which we have identified is Qir2ir° . Its 
branching ratios from the \P\p and 2P\$ states are approximately an order of magnitude 
larger than those for the 67r mode. Modes with charged pions and an odd number of neutral 
pions are forbidden by G-parity unless subsystems contain isospin-violating decays such as 
i] — > 7T + 7r~7r°. Indeed, 67T7T and 67r37r° decays are not seen at a statistically significant level. 
The 6tcAtt° mode involves fourteen particles, while we consider modes with a maximum of 
twelve. 

(2) The branching ratios for 87r27r° states from 1P\^ an d 2P 12 also exceed those for 8n 
by a considerable margin. Again, G-parity conservation explains why one does not see a 
significant signal for 8tt7t°. 

(3) Modes with one or more KK pairs in addition to charged pions are exempt from the 
G-parity selection rule because a KK pair can have either G-parity. 

(4) The 47r27r° mode has a larger significance than either 4n or AttAtt . Typically in the 
decay of an isospin-zero particle one should expect to see the same number of n + , n~ , and 
7T° [lOj, and this is reflected to some extent in individual modes. 

(5) The 14 modes constitute a total of less than a percent of all expected hadronic modes 
of the 1Pi,2 states. The ability to identify even such a small subset of the lPi^ hadronic 
decays depends to a large extent on CLEO's ability to reconstruct one or more neutral pions. 
Using only charged tracks one would reconstruct an order of magnitude fewer decays. 

We have also studied the El transitions T(3S) — > jXbO-Pj)- These transitions are sup- 
pressed by small overlaps of wave functions in the dipole matrix element (lP\f\3S) (ll| . We 
investigate the ratios 

_ g[T(3S)-> 7X6(1^)] m 
J B[T(2S) - lXb{lPj)] ' 1 ) 

where the ratios of branching fractions are determined from fitted yields in T(3S) and T(2S) 
data, respectively, corrected for small differences in signal efficiencies. In modes with any 
neutral pions or 7] mesons, substantial backgrounds arise from subsequent ^^(lPj) — > 7T(15) 
decays, whose photons are similar in energy to those in T(3<S') — > 7Xfc(lPj)- To eliminate 
such backgrounds, we restrict attention to 167 modes not involving tc° or rj mesons. Fig. 
|3]^a) shows the M inv distribution based on the sum of such modes in T(2S) data. The 
background is represented by the exponential of a polynomial, fitted to a shifted invariant 



6 



TABLE I: Reconstruction efficiencies for T(nS) — > 7 + Xb(( n — —> + X (e in units of 10 2 ), event yields (N), and signal significances 

(a) for each of the transitions for each of the 14 modes. In this and subsequent tables n = n , K = K ± . 
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TABLE II: Values of B[T(nS) -» 7x 6 ((ra - l)Pj)] x B[x b {{n - l)Pj) -> Xi] (10~ 5 ). Upper limits at 90% C.L. are set for modes with less 
than 3 a significance (see Table [J). 
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TABLE III: Values of B[xb(( n — ~^ Xi] (10 4 ). Upper limits at 90% C.L. are set for modes with less than 3 a significance (see TableU]). 
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< 10.8 


8.6 ±2.0 ±2.4 


9.6 ± 2 


.3 ±2.8 


3.9 ±1.6 ±1.1 


4.7 ± 1.8 ± 1. 


.4 




< 1.7 


< 6.7 


9.2 ±2.3 ±2.5 


6.7 ±1 


.9 ± 1.9 


< 5.0 


< 4.5 




6tt 


< 0.8 


< 0.7 


1.8 ±0.4 ±0.4 


1.2 ±0 


.3 ±0.3 


0.7 ±0.3 ±0.2 


0.9 ± 0.3 ±0. 


.2 


67T27T 


< 5.9 


< 12.3 


17.2 ±2.7 ±4.8 


11.9 ±2 


.4 ±3.4 


10.2 ±2.2 ±2.8 12.1 ±2.5 ±3. 


.6 


6tt2K 


2.4 ±0.9 ±0.7 


< 1.5 


2.6 ±0.6 ±0.7 


2.0 ±0 


.6 ±0.5 


< 0.8 


1.4 ± 0.5 ±0. 


.4 


6tt2K1tt° 


< 9.9 


< 7.3 


7.5 ± 1.6 ±2.1 


6.1 ± 1 


.4 ±1.8 


3.7 ±1.2 ±1.0 


4.2 ± 1.2 ±1. 


.2 


8vr 


< 0.7 


< 1.7 


2.7 ±0.6 ±0.7 


1.7 ±0 


.5 ±0.5 


0.8 ±0.4 ±0.2 


0.9 ± 0.4 ±0. 


.3 


87r27r° 


< 20.5 


< 6.5 


14.0 ±3.5 ±4.3 19.2 ±3 


.7 ±6.0 


18.5 ±4.4 ±5.6 12.6 ±3.5 ±4. 


.1 




FIG. 3: Comparison of spectra, based on 167 modes not involving 7r° or rj mesons, in transitions 
from (a) T(25) — > 7Xb(l.P) and (b) T(3S) — > -yxbO-P)- The background shape in (a) is derived 
from T(15) data as described in the text. The dashed lines in (a) are fitted Xb( n Pj) signals of 
each J state while those in (b) correspond to 90% C.L. upper limits on signal yields. 



mass distribution of T(IS') data (smooth dotted curve), and the peak positions are fixed at 
the known masses 

Using the XbO-Pj) signal shapes obtained in this fit, we fit T(3S) data with a constant 
background, as shown in Fig. [3](b). The fit gives statistical significances of 0.7cr, 0.0a, and 
2.6a for signals consistent with the transitions T(3S) — > 7X&(lPo,i,2), respectively. We find 
R 2 = (15.7 ±7.6 ±2.2) x 10~ 2 and 90% C.L. upper limits P(o,i,2)'< (21.9,2.5,27.1) x 10~ 2 . 
Using known values of B[T(25) -> 7X b (lP)] [1], we then find B[T(35) -> 7X b (lP 2 )] = 
[11 ± 6(stat.) ± 2(syst.) ± 1] x 10~ 3 , where the last uncertainty comes from B[T(2S) — > 
7X6(1-^2)] [l|- While most systematic uncertainties are canceled in the ratio of yields, our 
total uncertainty (14%) is dominated by variations in the signal as we change the width of 
the range over which we fit T(3S) decays. Our nominal fit range is 9800-9990 MeV, varied 
to 9800-9950 MeV and 9750-10050 MeV. Although this variation is well within statistical 
fluctuations, we conservatively take it as a possible systematic uncertainty. 

We set 90% C.L. upper limits £[T(3S) -> 7X&QLP0)] < 9.2 x 10~ 3 , consistent with the 
value of (3.0 ± 0.4 ± 1.0) x 10~ 3 reported in Ref. 3, B[T(35) -> !Xb(lPi)] < 1.9 x lO" 3 , 
and B[T(3S) 7X6(1P 2 )] < 20.3 x 10 -3 . Our results are compared with some theoretical 
predictions in Table ITVl 

We have presented the first observations of decays of Xb(lPj) and Xb(2Pj) to exclusive 
final states of light hadrons. These results can be of use in validating models for frag- 
mentation of heavy states, and in searching for states of mass ~ 10 GeV/c 2 via their 
exclusive decays. We also find upper limits for the rates of the suppressed El transitions 
T(3S0 - 7X & (lPo,i, 2 ). 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
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TABLE IV: Comparison of measurements and predictions 12|] for suppressed El transition rates 
in units of eV. Experimental measurements are based on r-w 3S ) = 20.32 keV [l|. 





J = 


J = 1 


J = 2 


Inclusive expt. [4] 


61 ±22 


- 


- 


Exclusive expt. (this work) 


< 186 


< 38 


< 413 


Moxhay-Rosner (1983) 


25 


25 


150 


Grotch et al. (1984) (a) 


114 


3.4 


194 


Grotch et al. (1984) (b) 


130 


0.3 


430 


Daghighian-Silverman (1987) 


16 


100 


650 


Fulcher (1990) 


10 


20 


30 


Lahde (2003) 


150 


110 


40 


Ebert et al. (2003) 


27 


67 


97 



(a) Confining potential is purely scalar. 

(b) Confining potential is purely vector. 



luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
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